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Abstract 

We present an analytical study of four-wave mixing instabilities in 
tapered fibers and photonic crystal fibers. Our approach avoids the use 
of Taylor expansion for the linear susceptibility and the slowly- varying 
envelope approximation. This allows us to describe the generation of 
sidebands strongly detuned from the pump wave with simultaneous ac- 
count for the entire dispersion characteristic of a fiber, which is found 
to be important for describing properly the key role of the parametric 
instabilities in the supercontinuum generation in these kind of fibers. 
OCIS code: 060.4370, 190.2620 

INTRODUCTION 

Four-wave mixing (FWM) is one of the most fundamental processes in optics, 
consisting in the generation of a pair of Stokes and anti-Stokes photons out 
of two strong pump photons [JJ, . FWM can be observed in a wide range of 
materials with Kerr nonlinearity including optical fibers 0. The efficiency 
of FWM (which is a coherent process) strongly depends on the so-called 
phase matching conditions imposed on the wave-vectors of the pump and 
generated waves ||. In a conventional single-cladding optical waveguide, 
the phase matching conditions are satisfied only in the vicinity of the zero- 
dispersion point, as it has been demonstrated for dispersion-shifted fibers 
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and for standard telecommunication fibers 0, [J, even if there are several 
techniques (for example using modal birefringence) which can extend the 
phase-matched frequency-range 0. 

Several experimental results on observation of four-wave mixing in pho- 
tonic crystal fibers (PCFs) @ and tapered fibers (TFs) [f| have been pub- 
lished over the last several years [R], 11]. These relatively new kind of 
fibers represent an ideal system for investigating the optical nonlinearities in 
fused-silica, because of their unique dispersive and nonlinear properties (see 



for example |L2], |13|)- in particular, the enhanced nonlinearity due to the 
smallness of the effective core area can increase dramatically every nonlinear 
effect. 

The phase matching conditions for these fibers are found to be quanti- 
tatively different with respect to ordinary fibers: phase matching in PCFs 
and TFs can be achieved for a long range of pump wavelengths, because the 
strong waveguide contribution to the overall dispersion permits a compen- 
sation of the material dispersion for a broad window of frequencies. On the 
other hand, the improved nonlinearity can generate a nonlinear coefficient 
which can further improve the compensation in the phase. 

Moreover, the dispersion characteristics of TFs || can be very similar 
to the ones found in PCFs. This is because in a typical index-guiding PCF 
most of the light is guided in the tiny silica core surrounded by a periodic 
structure of large air filled holes separated by thin membranes, which makes 
it very similar to TFs with guiding due to total internal reflection at the 
interface between silica and air |J, and in fact the eigenvalue equation used 
to find the group velocity dispersion (GVD) in TFs can describe with a good 
approximation the propagation of light in PCFs as well |TJ]], provided that we 
substitute the complex cladding structure of PCF with an effective refractive 



index [TS]. Thus, because of the fact that the GVD profile is qualitatively 



similar in these fibers, they can have analogous properties in the dynamical 
behavior of light propagation, as we shall see in the following. 

Fig. 1 shows numerically calculated examples of the frequency depen- 
dence of the dispersion parameter D = — (27ru 2 /c)d 2 j3, where (3 is the prop- 
agation constant of the fundamental mode 0], for a PCF with an elliptical 
core surrounded by 4 periods of air holes with an average diameter of 0.84 fim 
and a pitch of A = 1.46 fim, and a TF with 1 fim of core diameter, used in re- 
cent experiments on solitonic fission and super continuum generation (SCG) 



carried out by our group [0, ITB, I7J. Both dispersion characteristics show 



two zero dispersion (ZD) points limiting a central region of anomalous group 
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velocity dispersion (GVD), which is a typical situation for these fibers. 

In several important studies FWM in PCFs has been identified as one of 
the primary mechanisms contributing to generation of optical supercontin- 
uum by ultrashort pulse propagation, at least in the first stage of the spatial 
propagation g [H| pf. 

However, the analytical calculations complementing these results were 
either based on the standard analysis of some improved versions of the non- 
linear Schroedinger equations (NLSE) |19[ or on the direct analysis of the 
phase matching conditions || [L(| [H], |2(| . 

In this paper we investigate theoretically the role of parametric instabili- 
ties in PCFs and TFs, and we establish a deep link between these instabilities 
and the broad spectrum obtained in SCG, stimulated by recent experimental 



works (II |T0[) in which the primary mechanism of spectral broadening is 



identified as the combined action of stimulated Raman scattering and para- 
metric FWM. We shall find out that the unusual GVD profile of these kind 
of fibers will play a key role in determining the instability regions, and the 
analysis of these instabilities cannot be rigorously treated with a Taylor ex- 
pansion around the carrier frequency of the pulse because of its flatness and 
curvature. Moreover, we propose a mechanism for exciting a broad window 
of frequencies, based on the merging between different instability regions. 



MODEL EQUATION AND STABILITY ANALYSIS 

The formal mathematical approach to the FWM problem is to study the 
stability of the strong pump wave in the presence of weak perturbations with 
different frequencies. The typical instability of the pump wave occurring in 
optical fibers is the so-called modulational instability (MI) ||, which is a 
particular case of the more general phenomenon of parametric FWM. More 
specifically, MI can be interpreted as a four-photon process driven by the 
self-phase modulation (SPM), and it is a general feature of wave propagation 
in dispersive nonlinear media. In optical fibers MI is due to the conserva- 
tive interaction between nonlinearity and linear dispersion, which leads to 
a modulation of amplitude and phase of a continuous wave (CW) in pres- 
ence of quantum noise or any other weak perturbation. Traditionally the 
theoretical approach to describe MI in fibers is performed using the non- 
linear Schroedinger equation, that is the simplest equation which takes into 
account the interplay between dispersion and nonlinearity. NLSE in the con- 
text of fiber optics is derived from Maxwell's equations under the assumption 
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of the weakly-guiding weakly-nonlinear approximations and using a reduc- 
tive perturbation procedure consisting of the spatio-temporal slowly-varying 
envelope approximation (SVEA) and Taylor expansion of the linear suscepti- 
bility: x 0- The expansion of the latter is often done up to second order, 
which corresponds to the quadratic frequency dependence of the mode prop- 
agation constant (3: (3 ~ u 2 . If GVD at the carrier frequency u) is relatively 
small, then higher order terms in the expansion of x become relevant to 
describe propagation of waves in fibers in general and MI in particular, see 
e.g. 
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However, for some fibers with complex dispersion profile and/or in cases, 
when FWM results in generation of frequencies far detuned from the pump, 
the NLSE-based approach and in particular the Taylor expansion for x 
become inapplicable and new approaches have to be developed. One of them 
is to introduce a system of coupled NLSEs for pump, Stokes and anti-Stokes 
waves ||. 

Another more sophisticated method is to avoid the Taylor expansion of 
X^ and analyze the frequency mixing using a single equation, see p2[ and 



23 . In p2f this was mainly done in order to analyze the influence of the 



Raman effect on MI, and in it was applied to describe FWM in dispersion- 
flattened fibers. Our goal here is to further develop the theories reported 
22J, |2B| and consider PCFs and TFs as examples, thereby providing a 



m 



theoretical underpinning for recent experimental results || [T(], [TTJ, [L9| . 

We begin our analysis from the nonlinear wave equation derived directly 
from Maxwell's equations: 



1 

c 

where the linear polarization is defined as 



V 2 E - V(V • E) - -d 2 (E - P L - P NL ) = 0, (1) 



Pl = J x (1) (t-t>,x,y)E(t')dt>. (2) 

E is the electric field, and x^ is the linear susceptibility of the dielectric 
medium (in our case silica glass), which depends on time and the trans- 
verse coordinates x, y. Pnl is assumed to have a simple Kerr form: P NL = 
X3|E| 2 E. Even though the Raman effect (i.e., a delayed response of the non- 
linear susceptibility) can be incorporated into the theory developed below, 
analysis of details of its influence on FWM is well known and goes beyond 



the objectives of this paper [g^, [24[ . 
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Our next step is to reduce Eq. ([!]) to an equation in z and t only. In order 
to achieve this we first transform Eq. (|l|) from time to frequency domain, 
using Fourier transform J 7 , and then separate transverse and longitudinal 
degrees of freedom through the approximate factorization 

JFE(z, y, z, t) = E(x, y, z, u) ~ F(x, y, w)E{z, u). (3) 

F is the fundamental eigenmode of the linear waveguide having propagation 
constant (3{u). (3{uj) incorporates both material and waveguide contributions 
into the overall fiber dispersion. The factorization becomes possible due to 
the weakness of the energy transfer from the fundamental mode to the higher 



order modes 2^]. 



It can be shown that after using Eq. (|3|) the dynamics of the inverse 
Fourier transform T~ X E = E(z,t) of the amplitude E is governed by the 
equation 



JdlE - %E = dl 



f X { ejf(t-t')E(t')dt' + X 3\E\ 2 E , 



where x^ff ^ s an effective linear susceptibility of the fiber, with its Fourier 
transform given by 

= T X % = /3V/o; 2 - 1 (5) 

and 

JI F | 4 ^ 

* 3 = JWNs X3 ' (6) 



where S is the fiber area [H. %3, with good accuracy, can be considered as 
a frequency independent coefficient, which follows from the weak frequency 
dependence of F and X3 • 

Thus we have derived a version of the wave-equation with delayed linear 
and instantaneous nonlinear responses. Note that Eq. (f|) is an equation 
for the field, and not for the envelope. Expanding Xeff ^ n ^° Taylor series 
for small frequency detunings and making standard SVEA (in space and in 
time), one can easily reduce (§) to the conventional NLSE. 

We assume now that our fiber is pumped by a monochromatic continuous 
wave (CW) 

E = E e i{kz - Wt) + c.c, (7) 
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where E is a constant amplitude, k and u are the wave-number and the 
frequency of the pump wave. It is easy to see that Eq. (0) solves Eq. (f|) 
and represents a steady-state solution, provided that the following dispersion 
relation is satisfied: 



k 2 (u,\E \ 2 ) = ^ 



X^M + xsl^ol 2 • (8) 



To find the frequencies generated within the fiber as a result of the desta- 
bilization of the pump wave, we perturb solution (|7|) with small complex 
perturbations e(z,t): 

E =[E + e(z,t)}e l( > kz - UJt) + c.c. (9) 

After substitution of (§) into the governing equation (f|) we disregard all the 
terms nonlinear in e. Then we assume that the complex perturbation has 
the following simple form, which excites both the Stokes and the anti-Stokes 
bands: 

e(t) = Aexp(inz - i8t) + B*exp(-iK*z + iSt), (10) 

and derive a linear eigenvalue problem for the vector (A, B) T with an eigen- 
value k. It is possible to prove that the solvability condition for A and B 
requires that the following fourth-order polynomial equation in k is satisfied: 

- (fej + kl + 2k 2 )n 2 + 2k(k 2 + - k 2 _)n 

+ (kl ~ k 2 )(k 2 _ - k 2 ) - |U V + S) 2 (u - 5) 2 = 0, (11) 

where k± are the effective wave numbers of the Stokes and anti-Stokes waves, 
which include nonlinear corrections: k± = k(u ± (5, 2 1 _ZE7 1 2 ) - 

Solution (|?p becomes unstable provided that Eq. (]TT|) has at least one 
root such that Jtti(k) < 0. As usual, the instability gain is given by 

g(5)=2\Im(K)\. (12) 

In the next section we explore the numerical solutions of Eq. flTT|) . 

FOUR- WAVE MIXING INSTABILITIES 

We have scanned numerically the dependence of the imaginary parts of all 
four roots of (|TT|) from S for different values of pump frequency to and pump 
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power choosing the dispersion profiles shown in Fig. 1. Typical instability- 
gain profiles that we have obtained are shown in Fig. 2. We have found 
that there are two - symmetric with respect to the pump frequency - pairs 
of instability peaks in the spectrum of the perturbation. All these instability 
peaks are associated with the same root k = n cr of Eq. (JTTJ) . The regions 
nearest to the pump frequency can be traced back to the conventional MI, 
known from the NLSE, and the second peaks are due to the FWM process 
with far detuned frequencies, as we shall see below. 

Condition for both of these instabilities to exist is that the pump fre- 
quency should lay between the two ZD points or slightly outside this interval, 
see Figs. 2 and 3. 

In Fig. 2(a,c) the instability gain for our PCF and TF is shown, pumping 
in anomalous dispersion, with a frequency in between the two ZD points, 
ujzdi and ujzd2- The GVD profiles are also depicted [curve (1)], in order 
to show that the magnitude of the instability detunings is of the order of 
the anomalous regime width, \uzd2 ~ &zdi\- Moving slightly towards normal 
dispersion [Fig. 2(b,d)], a gap between the central frequency and MI appears. 
This has been observed in several numerical simulations, for example in Coen 
et al. and in Dudley et al. |H| . However, increasing the pump power one 



can observe a merging phenomenon between the two kinds of instabilities 
[curve (3) in Fig. 2], which is even more evident in normal dispersion, due to 
the fact that in this regime a slight change in the pump frequency towards the 
deep normal dispersion region corresponds to a large change in the position 
of MI, while the FWM peaks are almost fixed. 

There is a simple theoretical explanation for having instabilities even 
when ft" > 0, using the conventional phase matching condition 



k(u + 5, 2\E \ 2 ) + k(u - 6, 2\E \ 2 ) - 2k{u, \E \ 2 ) = 0. (13) 

Expanding k in Taylor series, knowing that ft (to) = (aye) 1 + X e Jf(u}) , 
and considering the detuning 5 as a small parameter, one obtains that there 



is a cancellation between the odd order derivatives in the expansion ||21|| , and 
it can be easily proved that Eq. ( |I3| ) can be written as 

±ft»6 2 +±ft""6 4 +... = U 2 f [ft"{u + £5) + ft"{u - £*)] (1-Odf = -xsl^ol 2 

«/ 

(14) 

Even if the pump is injected in a condition of small normal dispersion (ft" > 
0), there can be a contribution due to the even order terms present in Eq. 
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So in principle the integral in Eq. ( fH] ) can be negative, making possible 
the generation of instability in the normal dispersion regime. This simple 
analysis of course is valid for small frequency detunings 5, but our equation 
(|TT| ) does not have this limitation. 

One point we want to stress is that the dispersion inserted in Eq. (|TTD has 
to be considered globally, and not locally, in the sense that for a given pump 
frequency the entire shape of the GVD curve between and outside the two 
ZD points has significant influence on the FWM instabilities. Thus, Taylor 
expansion around a single ZD point is not adequate in this case if one wants 
to capture all the features of the process. 

Finding zeros of the function dsIm(K cr ), which correspond to the positions 
of the maximal gains occurring for 5 = 5 cr , we calculate a diagram showing 
frequencies ujas,s = uj±5 cr of the generated waves as functions of u, see Fig. 
3. Fig. 3(a) shows the position of Stokes and anti-Stokes waves for every 
pump frequency and for three different powers in the case of PCF, while 
Fig. 3(b) is for our TF. Dashed curves correspond to the two ZD points. 
It is possible to see the merging process when increasing the pump power 
(respect to a fixed pump frequency) or moving the pump frequency (respect 
to fixed pump power). 

Note that increasing the power the distance between the MI peaks in- 
creases considerably, while the position of FWM peaks remains approxi- 
mately the same. Moreover, one of the two ZD points (in our case the 
blue-shifted one, but of course it depends on the particular GVD that we 
choose) is more "efficient", in the sense that the frequency range excited by 
the instabilities is larger respect to the other ZD point. An important dif- 
ference (due to the absence of a Taylor expansion in our equations) between 
Fig. 3 and other figures present in the literature (for example Ref. |9|, [1Q|] ) 
is that our curve is closed, and then the unstable region in the normal dis- 
persion is finite. A simple prediction which follows from this consideration 
is that it should exist a pump frequency in the normal regime for which all 
these instabilities are strongly suppressed and completely disappear. 

To demonstrate that instabilities found here are essentially FWM insta- 
bilities we also show the frequencies of the generated waves calculated using 
the standard phase matching condition flT3"|). One can see remarkable agree- 
ment between results obtained using the exact solution of Eq. (|ll|) and 
condition (0), confirming the fact that starting with a forward pump wave, 
we can neglect the influence of the back-propagating wave in the perturba- 
tion analysis. Note that the phase matching diagrams presented in |], 1(| 
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in the vicinity of the ZD point with largest frequency capture only a pair of 
instability peaks associated with conventional MI. 

Existence of the far detuned FWM branches shows that contribution of 
the FWM processes to the generation of optical super continuum can be in 
some cases even stronger than it has been so far anticipated. Also the initial 
stage of super continuum generation can be attributed to the threshold mech- 
anism of merging between the FWM peaks and the MI peaks [|y|, clearly 
visible from Fig. 2 when one increases the pump power. Through this mech- 
anism it is possible to excite all the frequencies between the two FWM peaks, 
and it may represent an important piece towards the understanding of SCG. 

Making qualitative comparison of our results with these experiments is 
however, complicated by the silica absorption, which is not taken into ac- 
count in our theory, that is substantial in the far infrared (over 2 microns). 
Therefore one should expect that far red detuned radiation is strongly ab- 
sorbed, and this can potentially modify instability conditions and shift the 
position of the frequencies with maximal gain. Even the Raman effect can 
affect the instability gain introducing an asymmetry, i.e. a suppression of 
blue-shifted frequencies respect to the red-shifted ones. 

Approximating the propagation constant in Taylor series, one can show 
that secondary FWM peaks appear in correspondence of the frequency de- 
tuning 

-12/3" (u)^ 1/2 



m - h^f J (15) 

From Eq. fllS]) it is evident that the second-order dispersion /3" and the 
fourth-order dispersion f3"" must have different sign in order to generate the 
peaks associated to FWM. 

Thus pumping the fiber in the region where GVD as function of u is 
flatter, i.e. f3""{u) is smaller, will generate a broader spectrum. In our 
examples f3""{uo) is smaller near the blue shifted ZD point, therefore pumping 
in its vicinity one can expect generation of the broadest possible spectrum. 

An important estimate, which can be found analytically, is for the critical 
power required for merging of the MI and FWM peaks. It can be shown that 
the power required for merging is given by 

2 ^ 6n e// Mc|^M| 2 

m ~ Xsuj |/3""(u;)| l ° j 

is less for pump frequencies close to the ZD points. Moreover, because 
of the factor 1/cu in this relation, E^ is less for the blue-shifted ZD point 
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respect to the red-shifted one. It follows that in order to improve the width 
of the generated spectrum one needs to find a compromise between the power 
required to activate the process, which depends on |/3"| 2 /|/3""|, and a small 
value of | at the pump frequency, which regulates the relative distance 
between the peaks involved in the mechanism. Thus flat and small GVD 
are the two factors responsible for generation of broad FWM spectrum at 
low powers, which agrees with experimental finding on generation of optical 
supercontinuum [g g |H], 0, [Hj, [19], £§• 

The plots of the FWM gain in Fig. 2 are superimposed on the plots 
showing the dispersion profile of the fiber. It is clear from this figure that 
dispersion varies quite dramatically through the relevant frequency domain; 
moreover, the order of magnitude of the detuning relative to the generated 
frequencies is comparable with the pump frequency itself, such that \S\/u> can 
not be considered as a small parameter. These are the two reasons which 
validates the necessity of avoiding Taylor expansion for x ■ 

CONCLUSIONS 

Parametric instabilities in PCFs and TFs have been studied in the context 
of a general wave equation, avoiding the usual approximations in optics, like 
the Taylor expansion of the linear susceptibility and the suppression of the 
contribution of backward waves. Other regions of instability has been found, 
corresponding to degenerate FWM, which can exist even in a small region 
of normal dispersion. These peaks of instability are found to be relatively 
close to the pump frequency in our examples of PCF and TF, but still too 
far to be detected through ordinary methods of theoretical and numerical 
analysis. We have evaluated analytically the power required to merge this 
regions to the modulational instability regions, and we believe that this mech- 
anism (regulated by the usual fourth order dispersion (3"" in the language of 
the ordinary NLSE) can be of fundamental importance in the description 
of supercontinuum generation in fibers which have a GVD with two zero 
dispersion points. 

ACKNOWLEDGMENTS 

We acknowledge J.K. Knight, P.St. J. Russell and T.A. Birks for several useful 
discussions. 



10 



References 

[1] R.H. Stolen, and G.D. Bjorkolm, IEEE J. Quantum Electron. 18, 1062 
(1982). 

[2] O. Aso, M. Tadakuma, and S. Namiki, Furukawa Rev. 19, 63 (2000). 

[3] G.P. Agrawal, Nonlinear Fiber Optics (Academic Press, San Diego, 
2001). 

[4] W. Washio, K. Inoue, and S. Kishida, Electron. Lett. 16, 658 (1980). 

[5] D.K. Serkland and P. Kumar, Opt. Lett. 24, 92 (1999). 

[6] R.H. Stolen, M.A. Bosch, and C. Lin, Opt. Lett. 6, 213 (1981). 

[7] J.C. Knight, T.A. Birks, P.St. J. Russell, and D.M. Atkin, Opt. Lett. 21, 
1547 (1996). 

[8] T.A. Birks, W.J. Wadsworth, and P.St. J. Russell, Opt. Lett. 25, 1415 
(2000). 

[9] S. Coen, A.H.L. Chau, R. Leonhardt, J.D. Harvey, J.C. Knight, W.J. 
Wadsworth, and P.St. J. Russell, J. Opt. Soc. Am. B 19, 753 (2002). 

[10] J.M. Dudley, L. Provino, N. Grossard, H. Maillotte, R.S. Windeler, B.J. 
Eggleton, and S. Coen, J. Opt. Soc. Am. B 19, 765 (2002). 

[11] J.E. Sharping, M. Fiorentino, A. Coker, P. Kumar, and R.S. Windeler, 
Opt. Lett. 26, 1048 (2001). 

[12] A.M. Zheltikov, Physics Uspekhi 43, 1125 (2000). 

[13] J. Broeng, D. Mogilevstev, S.E. Barkou, and A. Bjarklev, Opt. Fiber 
Tech. 5, 305 (1999). 

[14] A.W. Snyder, and J.D. Love, Optical Waveguide Theory (Chapman and 
Hall, London, 1983). 

[15] J.C. Knight, T.A. Birks, P.St. J. Russell, and J. P. de Sandro, J. Opt. 
Soc. Am. A, 15, 748 (1998). 



11 



[16] A. Ortigosa-Blanch, J.C. Knight, and P.St. J. Russell, J. Opt. Soc. Am. 
B (2002, to be published). 

[17] W.J. Wadsworth, A. Ortigosa-Blanch, J.C. Knight, T.A. Birks, T.P.M. 
Man, and P.St. J. Russell, J. Opt. Soc. Am. B, 19, 2148 (2002). 

[18] A.V. Husakou, and J. Herrmann, J. Opt. Soc. Am. B 19, 2171 (2002). 

[19] G. Millot, A. Sauter, J.M. Dudley, Opt. Lett. 27, 695 (2002). 

[20] S. Song, C.T. Allen, K.R. Demarest, and R. Hui, IEEE J. Light. Tech. 
17, 2285 (1999). 

[21] M.J. Potasek, Opt. Lett. 12, 921 (1987). 

[22] K.J. Blow, D. Wood, IEEE J. Quantum Electron. 25, 2665 (1989). 

[23] M. Yu, C.J. McKinstrie, and CP. Agrawal, Phys. Rev. E 52, 1072. 
(1995). 

[24] E.A. Golovchenko, and A.N. Pilipetskii, J. Opt. Soc. Am. B 11, 92 
(1994). 

[25] A.V. Husakou, and J. Herrmann, Phys. Rev. Lett. 87, 203901 (2001). 

[26] J. Herrmann, U. Griebner, N. Zhavoronkov, A. Husakou, D. Nickel, J.C. 
Knight, W.J. Wadsworth, P.St. J. Russell, and G. Korn, Phys. Rev. Lett. 
88, 173901 (2002). 



12 



List of figures 




Figure 1: Group- velocity dispersion (GVD) for our examples of PCF (a) and 
TF (b) described in the text. 
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Figure 2: Curves (1) show the dependence of the rescaled parameter D(uS) 
on frequency for the chosen PCF and TF with 1 fim core diameter. To get 
physical values of D one needs to multiply the values for the instability gain 
by the factor 17 ps/nm/km. Curves (2) and (3) show the dependence of the 
instability gain Im(hz cr ) for PCF (a,b) and TF (c,d) on frequency for different 
pump powers P = \E \ 2 . (2) is for \E \ 2 < \E m \ 2 and (3) is for \E \ 2 > \E m \ 2 . 
The central frequency in every picture corresponds to the pump frequency. 
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Figure 3: Diagrams showing positions of the maximal instability gain vs 
pump frequency calculated from Eq. (7) for different powers P = \E \ 2 in 
PCF (a) and TF (b). Dots show the same values calculated using phase 
matching condition (8). ZD points are indicated by dashed vertical lines. 
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